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The 2 H(e,e'p)n cross section at a momentum transfer of 3.5 (GeV/c) 2 was measured over a 
kinematical range that made it possible to study this reaction for a set of fixed missing momenta as 
a function of the neutron recoil angle nq and to extract missing momentum distributions for fixed 
values of nq up to 0.55 GeV/c. In the region of 35° < nq < 45° recent calculations, which predict 
that final state interactions are small, agree reasonably well with the experimental data. Therefore 
these experimental reduced cross sections provide direct access to the high momentum component 
of the deuteron momentum distribution in exclusive deuteron electro-disintegration. 
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The understanding of the short-range structure of the 
deuteron is of fundamental importance for the advance- 
ment of our understanding of nuclear matter at small 



* deceased 

t present address: Indiana University, Bloomington, Indiana 47405 
tpresent address: Renaissance Technologies LLC, East Setaukct, 
New York 11733 

^present address: Massachusetts Institute ol Technology, Cam- 
bridge, Massachusetts 02139 



distances. To probe the short-range properties of the 
deuteron, one has to investigate configurations where the 
two nucleons come very close together and are strongly 
overlapping. The basic problem is to what extent these 
configurations can be described simply in terms of two 
nucleons with high initial relative momenta. The ulti- 
mate quantity to be investigated in this case is the high 
momentum component of the deuteron wave function. 
Traditionally three classes of reactions are used to study 
the high momentum part of the deuteron wave func- 
tion: elastic scattering, inclusive and exclusive electro- 
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disintegration reactions. 

Elastic electron-deuteron scattering probes the in- 
tegrated characteristics of the wave function via the 
deuteron form-factors. At large 4-momentum trans- 
fer, — Q 2 , the scattering becomes sensitive to the high 
momentum component of the deuteron wave function. 
The analysis of experimental data pQ showed that, at 
presently available energies, it is practically impossi- 
ble to discriminate between different theoretical ap- 
proaches 0[3] used to calculate the deuteron elastic form- 
factor A(Q 2 ). One needs additional constraints on the 
deuteron wave function at short distances. 

Inclusive, quasi-clastic (e, e') reactions provide an- 
other way of probing high momentum components of 
the deuteron especially at high Q 2 and in the xb > 1 
region [JHZ] where xg — Q 2 /2Mv (M is the nucleon 
mass and v is the energy of the virtual photon) is the 
Bjorken scaling variable. In this regime the cross section 
depends on an integral of the deuteron momentum dis- 
tribution with the longitudinal nucleon momentum com- 
ponent (with respect to the virtual photon momentum 
q) as the lower limit. However, the difficulties of ensur- 
ing small contributions from inelastic processes (growing 
with Q 2 ) and final state interactions (FSI) at large xb 
(see e.g. [HI IS]), reduce the sensitivity to the deuteron 
wave function at short inter-nucleon distances, although 
the high-momentum component is certainly important in 
this kinematics. 

The most direct way of studying high nucleon mo- 
menta is to investigate the quasi-elastic (QE) electro- 
disintegration of the deuteron via the 2 H(e, e'p)n reac- 
tion at high missing momenta (the momentum of the 
recoiling neutron) p m = q — pf, where pf is the mo- 
mentum of the outgoing, observed proton. Within the 
Plane Wave Impulse Approximation (PWIA) — p m cor- 
responds to the initial momentum of the target nucleon 
before the interaction. Thus the strategy in these studies 
is to probe the cross section at p m values as large as pos- 
sible. However, depending on the selected kinematics, 
these studies can be overwhelmed by final state inter- 
actions (FSI) where the outgoing proton interacts with 
the recoiling neutron, or by processes where the virtual 
photon couples to the exchanged meson (MEC) or where 
the nucleon is excited to an intermediate A state (IC). 
The dominance of FSI, MEC and IC has seriously af- 
fected previous experiments at Q 2 < 1 (GeV/c) 2 [TOtiM] 
leading to the overall conclusion that these experiments 
do not provide good constraints on the high momentum 
components of the deuteron wave function. 

The condition Q 2 > 1 (GeV/c) 2 is necessary in order 
to enhance contributions of reaction mechanisms which 
probe the short-range structure of the deuteron and to 
suppress competing long range processes for the follow- 
ing reasons: (i) the MEC contribution should be sup- 
pressed by an additional factor of (1 + Q 2 /A)^ 4 with 
A = 0.8 — 1 (GeV/c) 2 as compared to the QE contri- 
bution [T21 US] ; (ii) the large Q 2 limit should allow one 
to probe the wave function in the x > 1 region which 



is far from the inelastic threshold, thereby suppressing 
IC contributions; (iii) final state interactions of the out- 
going nucleon should follow the eikonal dynamics with 
a strong angular anisotropy dominating mainly in trans- 
verse directions. This situation generated a multitude of 
theoretical studies of the 2 H(e, e'p)n reaction in the high 
Q 2 regime [T51 [T7H2T] . The PWIA results of calculations 
described in ref. [TJl HH [H [251 [27] differ at larger p m due 
to differences in the wave functions and in details of the 
off-shell eletron nucleon interaction used, but all predict 
small FSI contributions (10 - 20%) for 35° < 9 nq < 45°. 

We report new 2 H(e, e'p)n cross sections measured at 
high momentum transfer for well defined kinematic set- 
tings. The wide range in missing momenta and neutron 
recoil angles allows one for the first time to access the 
high momentum components of the deuteron momentum 
distribution and probe the validity of current models of 
the reaction dynamics. The kinematic region covered 
overlaps with a recent 2 H(e, e'p)n experiment performed 
using CLAS at Jefferson Lab [25J which concluded that 
FSI and IC are dominating the momentum distribution 
except for p m < 0.1 GeV/c or 9 nq > 110°. However, 
to obtain reasonable statistical precision, the data were 
integrated over the full 9 nq range for the momentum dis- 
tributions and over a large p m -range for the angular dis- 
tributions in contrast to the data presented below. 

At a fixed Q 2 = 3.5 (GeV/c) 2 , the 2 H(e,e'p)n cross 
section was measured for specific missing momenta p m — 
0.1, 0.2, 0.4, 0.5 GeV/c, while the angle 6 nq of the re- 
coiling neutron with respect to q, was varied. 6 nq is also 
referred to hereafter as the recoil angle. For p m = 0.4, 0.5 
GeV/c, the largest recoil angles accessible were limited 
by the maximum momentum that the proton spectrom- 
eter was able to detect (3.1 GeV/c). Keeping p m and Q 2 
constant, required the energy transfer, the electron scat- 
tering angle, the proton final momentum and the pro- 
ton direction to be adjusted accordingly for each value of 
9 nq . As the energy transfer and recoil angle changed, Xb 
changed as well between 0.78 and 1.52, large xb values 
corresponding to small recoil angles. 

The experiment was carried out using the two high- 
resolution spectrometers in Hall A at Jefferson Lab at 
an electron beam momentum of 5.008 GeV/c. The left 
arm detected the electrons and the right arm the ejected 
protons. The deuterium target consisted of a 15 cm long 
cylinder filled with liquid deuterium and was part of the 
Hall A cryogenic target system [25]. An identical target 
cell filled with liquid hydrogen was used for calibration 
and to determine the coincidence efficiency. The electron 
beam was rastered over an area of 2 x 2 mm 2 and the 
liquids were continuously circulated in order to minimize 
density variations due to boiling. We found a typical 
reduction of the effective deuteron target thickness due 
to boiling by a factor of 0.94 ± 0.02 for an average cur- 
rent of 100 (x A. The cross sections were corrected for 
detector inefficiencies on a run-by-run basis and for an 
overall coincidence efficiency, determined from the mea- 
sured 1 H(e,e') elastic cross section which was found to 



3 



be 96.4 ± 2 percent of the published value from the fit of 
Table I in ref. [30 . The systematic error due to uncer- 
tainties in the measured kinematic variables were calcu- 
lated for each data bin and added quadratically to the 
statistical error. An overall error of 4.5 % was added 
to take into account errors in beam energy, beam charge 
measurements, detector efficiencies, target thickness and 
target boiling corrections. 

The spectrometer detection systems in the two arms 
were very similar: vertical drift chambers (VDC) were 
used for tracking and two scintillation counter planes 
(S1/S2) following the VDCs provided timing and trigger 
signals. In addition, the electron arm was equipped with 
a gas Cerenkov detector for electron/7r~ discrimination. 
We found that the gas Cerenkov detector was sufficient 
for the ir~ rejection in this experiment. At this large 
momentum transfer and at the large xb kinematics the 
7r~ background was not a concern. A detailed descrip- 
tion of the spectrometer systems and the target system 
can be found in reference |29j . The momentum accep- 
tance used for both spectrometers was set by software to 
S = Ap/po =±4 %, where po is the central momentum 
of the spectrometer. 

The solid angle of each spectrometer was defined by 
software cuts at the entrance of the first quadrupole 
magnet. In addition a second, global cut on the multi- 
dimensional acceptance of each spectrometer was ap- 
plied by means of R- functions [31 . The phase space 
acceptance was calculated using the Hall A Monte- 
Carlo code MCEEP [35]. The extracted cross sections 
were radiatively corrected using the Monte-Carlo code 
SIMC [331 GH] where the yield was estimated with a the- 
oretical calculation by J.M. Laget [55] that included final 
state interactions and reproduced the experimental yield 
quite well. 

Several spectrometer settings contributed to a full 
angular distribution and the condition of constant Q 2 
and constant p m was maintained for the central setting 
only. Within the phase space acceptance defined by the 
two spectrometers the kinematic variables varied slightly 
around their central values as a function of the angle 9 nq . 
This led to variations of the experimental cross section as 
a function of 9 nq that were independent of reaction mech- 
anism effects. In order to reduce these variations the ex- 
perimental cross section was divided by the PWIA cross 
section (JpwiA — n p(Pm)k<J cc i where k is a kinematic 
factor, np(p m ) refers to the Paris momentum distribution 
and a cc i is the de Forest CC1 off-shell cross section [35] 
calculated using the form factor parameterization from 
Table I of ref. [3D] . Theoretical cross sections were calcu- 
lated using the averaged kinematics determined for each 
bin in 9 nq , included bin center corrections and were di- 
vided by the same PWIA cross section as the experimen- 
tal ones. The resulting ratios R{9 nq ) — a exp j '<J pwia were 
then averaged for overlapping 9 nq bins and the resulting 
angular distributions are shown in Fig. [l] 

The angular distribution shown in Fig. [T^, for missing 
momentum p m = 0.2 ± 0.02 GeV/c shows a clear reduc- 
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FIG. 1: (Color online) The ratio R(0 nq ) = o- exp /a pw ia- 
(a) Pm = 0.2 GeV/c, (b) p m = 0.4 GeV/c and (c) p m = 
0.5 GeV/c. Solid (purple) lines (i) MS [H QS [27] us- 
ing the CD-Bonn potential, short dashed (green) lines (ii) 
JML [25J, dashed-double-dot lines JML with MEC and IC 
and long-dashed (orange) lines (in) JVO [22]. Insets: Rfsi = 
<W4t lo for 35° <6>„ 9 <45°. 



tion of R for 9 nq around 75° (xs~l) . For missing momen- 
tum p m = 0.4±0.02 GeV/c (Fig.^) andp m = 0.5±0.02 
GeV/c (Fig. [TJ;), R shows a peak at around 75° with a 
maximal value of ~ 1.6 and ~ 2.5, respectively. The de- 
pendence of R on 9 nq reflects the angular dependence of 
final state interactions at high Q 2 . At high energies, FSI 
are described in the eikonal regime where the fast proton 
rescatters off the spectator neutron which in turn recoils 
almost perpendicularly to q. 

This measurement confirms the strong angular varia- 
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tion of the 2 H(e, e'p)n cross section as a function of 9 nq 
observed in the CLAS experiment |28j for missing mo- 
menta 0.4 < p m < 0.6 GeV/c. It is the basic feature un- 
der consideration in the study of the color transparency 
phenomenon in few-body systems |36j . 

Three different theoretical calculations were obtained: 

(1) a calculation by M. Sargsian [T51 HH1 127] , referred to 
as MS below, using the CD-Bonn or the Paris potentials, 

(2) cross sections from J-M. Laget's model [25], referred 
to as JML, using the Paris potential and (3) results from 
the model of Jeschonnek and Van Orden 22 which will 
be labelled JVO below. The rclativistically covariant cal- 
culation of JVO is currently limited to small recoil angles. 

The model predictions of R = u C aic/ '&pwia are com- 
pared to experimental data in Fig. [T] For small missing 
momenta (p m < 0.2 GeV/c) and 9 nq < 30° all calcula- 
tions agree with each other within 20%. For larger an- 
gles and especially larger missing momenta, deviations 
between the experiment and the calculations and be- 
tween the different calculations themselves are consid- 
erably larger. For p m — 0.4 and 0.5 GeV/c MS correctly 
describes the rise of R with 9 nq , but predicts a faster fall- 
off after the maximum than is observed. JML predicts 
a considerably wider re-scattering peak. Including MEC 
and A excitation improves the agreement considerably 
at p m = 0.4, 0.5 GeV/c (Fig. flfc) but worsens the agree- 
ment for p m = 0.2 GeV/c (Fig. [TJi). The value of the 
maximum agrees with experiment for both calculations. 
For p rn = 0.4, 0.5 GeV/c and 9 nq < 50° MS and JVO de- 
scribe the data considerably better than JML excluding 
MEC and A contributions. The ratio Rf S i = o/si/CT*^ ia , 
where <J j s % is a calculated cross section including FSI and 
a pwia i s the corresponding PWIA cross section, demon- 
strates that all calculations show relatively small contri- 
butions of FSI for 35° < 9 nq < 45° (insets in Fig. [I}. For 
0.2 < p m < 0.5 GeV/c on average MS predicts between 
-18% and -5% FSI, JM between -9% and +5% and JVO 
between -12% and -16% FSI. This indicates a kinematic 
region where the cross section is dominated by PWIA and 
should therefore reflect directly the deuteron momentum 
distribution. 

We extracted the 2 H(e, e'p)n cross section for three 
sets of fixed angles 9 nq with a bin width of ±5°. For 
each of these three recoil angles, we determined the cross 
section as a function of missing momentum and calcu- 
lated the reduced cross section a re d = <J exp / {ka cc \) using 
the same form factor as previously (Table I in ref. [30] ) . 
We included all spectrometer settings that contributed 
to the same {p m ,9 nq ) bin and determined the averaged, 
reduced cross section as well as the averaged kinematics. 
The theoretical cross sections have been treated the same 
way as the experimental ones. The resulting experimen- 
tal momentum distributions are shown in Fig. [2j 

The results for 9 nq = 75° (Fig. ^) show a typical be- 
havior of reduced 2 H(e,e'p)n deuteron cross sections as 
a function of missing momentum, namely a 'flattening' 
of the cross section around p m > 0.3 GeV/c. This flat- 
tening has been observed in most previous measurements 



of the 2 H(e, e'p)n cross section at lower Q 2 which have 
been taken either at [13] xb ~ 1, where FSI dominates 
or at values xb < 1, where MEC and IC dominated [llj . 
From the measured angular distributions reported here 
and previously [2"8"] . we found that at this angle FSI con- 
tributions to the reaction are maximal. 

The experimental reduced cross sections are compared 
to a calculation by M. Sargsian with wave functions from 
the CD-Bonn (i) and the Paris (ii) potentials. The PWIA 
results are shown as solid curves and the ones includ- 
ing FSI as dashed (CD-Bonn) and short dashed (Paris) 
lines. Both calculations at 9 nq — 75° including FSI agree 
quite well with the measurement. The PWIA calcula- 
tions cannot reproduce the data for p m > 0.1 GeV/c and 
for p m > 0.45 GeV/c the two calculations increasingly 
deviate from each other. 

In contrast to 9 nq — 75° , the reduced cross sections 
9 nq = 35° and 9 nq = 45° , display a qualitatively different 
behavior as a function of p m . The fall off is considerably 
steeper for p m > 0.3 GeV/c and follows closely the gen- 
eral shape of the deuteron wave function in momentum 
space. At small 9 nq the calculated cross sections with FSI 
differ much less from the PWIA results and are sensitive 
to the type of NN potential used for p m > 0.45 GeV/c. 

We measured 2 H(e, e'p)n cross sections at a momen- 
tum transfer of 3.5 (GcV/c) 2 over a kinematical range 
that allowed us to study this reaction for a set of fixed 
missing momenta as a function of the neutron recoil angle 
9 nq . We experimentally confirmed the validity of the gen- 
eralized eikonal approximation which predicted a strong 
angular dependence of FSI contributions and kinematic 
regions where FSI contributions are small. 

The small kinematic bin size made it possible for the 
first time to determine missing momentum distributions 
for several, fixed values of the neutron recoil angle, 9 nq 
and to observe a qualitative change in their shape. With 
decreasing 9 nq the momentum distributions change from 
the typical form found in previous experiments to a shape 
that follows more closely the trend of the deuteron wave 
function in momentum space. This transition is consis- 
tent with decreasing FSI contributions and gives us for 
the first time a direct access to the high momentum com- 
ponent of the deuteron momentum distribution. We find 
that within the MS model the calculations using the CD- 
Bonn potential are in best agreement with the data. 

We acknowledge the outstanding efforts of the staff of 
the Accelerator and Physics Divisions at Jefferson Lab 
who made this experiment possible. This work was sup- 
ported in part by the Department of Energy under con- 
tracts DE-FG02-99ER41065, DE-AC02-06CH11357, the 
Italian Istituto Nazionale di Fisica Nucleare, the French 
Centre National de la Recherche Scientifique and the Na- 
tional Science Foundation. Jefferson Science Associates 
(JSA) operates the Thomas Jefferson National Acceler- 
ator Facility for the U.S. Department of Energy under 
contract DE-AC05-84ER40150. 
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FIG. 2: (Color online) The reduced cross section <T re< j(pm) as a function of missing momentum p m is shown in panels a, b and 
c for 9 nq = 35°, 45° and 75°, respectively, and a bin width of ±5°. CD-Bonn potential: dashed (blue) lines PWIA , solid (blue) 
lines FSI. Paris potential: dash-dot (green) lines PWIA, long-dashed (green) lines FSI. The PWIA results are for all angles 
identical. All calculations are from the MS model [TBI 1191 |27] . 
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